Soft and Supersoft X-ray Sources in Symbiotic 

Stars 1 



Chunhua Zhu a ' 6 , Guoliang Lii 6 , Zhaojun Wang a ' 5 and 

Jim Zhang * & 

a School of Science, Xi'an Jiaotong University, Xi'an, 710049, China 
b School of Physics, Xinjiang University, Urumqi, 830046, China 



Abstract 

Assuming that soft X-ray sources in symbiotic stars result from strong thermonu- 
clear runaways, and supersoft X-ray sources from weak thermonuclear runaways 
or steady hydrogen burning symbiotic stars, we investigate the Galactic soft and 
supersoft X-ray sources in symbiotic stars by means of population synthesis. The 
Galactic occurrence rates of soft X-ray sources and supersoft X-ray sources are from 
~ 2 to 20 yr _1 , and ~ 2 to 17 yr" 1 , respectively. The numbers of X-ray sources 
in symbiotic stars range from 2390 to 6120. We simulate the distribution of X-ray 
sources over orbital periods, masses and mass-accretion rates of white dwarfs. The 
agreement with observations is reasonable. 
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1 Introduction 



Symbiotic stars (SySs) are usually interacting binaries, composed of a cool 
star, a hot component and a nebula. The cool component is a red giant which 
is a first giant branch (FGB) or an asymptotic giant branch (AGB) star. The 
hot component is a white dwarf (WD), a subdwarf, an accreting low-mass 
main-sequence star, or a neutron star (Kenyon & Webbink, 1984; Miirset et 
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al., 1991). The peculiar optical spectrum of SySs is a very important and 
interesting phenomenon, and offers some exciting observational facts (Kenyon, 
1986; Miirset & Schmid, 1999; Belczyhski et al, 2000; Mikolajewska, 2003). 
Miirset et al. (1997) showed that the majority of the known galactic SySs are 
detectable X-ray sources. According to ROSAT observations of SySs, Miirset 
et al. (1997) divided X-ray sources into three distinct classes: a-type X-ray 
sources in which the photon energies are below 0.4 keV; /3-type X-ray sources 
in which the photon energies are roughly between 0.1 and 10 keV and the 
peak at about 0.8 keV; and 7-type relatively hard X-ray sources like GX 1+4. 
Following van den Heuvel et al. (1992), Yungelson et al. (1996) and Miirset 
et al. (1997), we call a-type X-ray sources as super-soft X-ray sources (SSSs), 
/3-type X-ray sources as soft X-sources (SSs) (Miirset et al., 1997). Due to hard 
X-ray sources resulting from accreting neutron stars (van den Heuvel et al., 
1992; Miirset et al., 1997; Orio et al., 2007), we do not discuss them in the 
present paper. 

In general, SSSs result from the steady hydrogen burning on accreting WD 
in binary systems. While for SSs, it was widely believed that they are arisen 
from the colliding winds from the giant and WD (Leahy & Volk, 1995; Miirset 
et al., 1997). In recent years, many theoretical studies on SySs have been 
published, e.g. Yungelson et al. (1995); Han et al. (1995a); Iben & Tutukov 
(1996); Hurley et al. (2002); Lti et al. (2006, 2007, 2008). However, there are 
few detailed theoretical studies about the Galactic population of X-ray sources 
in SySs. Yungelson et al. (1996) investigated the Galactic binary SSSs with 
white dwarf accretors by means of a population synthesis. They considered 
that SSSs in SySs result from the steady hydrogen burning and thermonuclear 
runaways, but they did not distinguish SSs from SSSs. 

In the present paper we model possible formation paths for the Galactic pop- 
ulation of SSs and SSSs in SySs with white dwarf accretors. In Section 2, we 
present our assumptions and describe some details of the algorithm. In Section 
3, we discuss the main results. In Section 4, the main conclusions are given. 



2 The model 



For binary evolution, we use a rapid binary star evolution (BSE) code of 
Hurley et al. (2002). Below we describe our algorithm from several aspects. 
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2.1 Symbiotic stars 



SySs are complex binary systems, and can be divided into two subgroups: 
'ordinary' SySs, which are assumed to burn hydrogen steadily; and symbiotic 
novae, which experience thermonuclear runaways in their surface hydrogen 
layers (Tutukov & Yungelson , 1976; Mikolajewska, 2003). There are many 
uncertain physical parameters which can affect the population of SySs. Lii et 
al. (2006) showed that the numbers of SySs and the occurrences of symbiotic 
novae are greatly affected by the algorithm of common envelope evolution, the 
terminal velocity of stellar wind v(oo) and the critical ignition mass AM^J 3 
which is necessary mass accreted by WD for a thermonuclear runaway. The 
structure factor of the stellar wind velocity aw and an optically thick wind 
give a small uncertainty. In this work, we use the models of SySs in Lii et al. 
(2006) and discuss the Galactic population of SSs and SSSs in SySs. 

Common envelope: Following Nelemans & Tout (2005) and Lii et al. (2006), 
we use two algorithms (a-algorithm and 7-algorithm) for common envelope 
evolution. The a-algorithm results in much shorter binary separation after un- 
dergoing common envelope phase than the 7-algorithm. In this work, we take 
the 'combined' parameter a ce A ce as 0.5 for a-algorithm. A ce is a structure pa- 
rameter that depends on the evolutionary stage of the donor. For 7-algorithm, 
7 = 1.75. 

v(oo): It is difficult to determine the terminal velocity of stellar wind v (00). 
Lii et al. (2006) used two completely different v(oo)s: 

(i) u(oo) = |w C sc, where t> e sc is the escape velocity. With ascent of star along 
giant branch, the stellar radius becomes larger and larger, which results in the 
lower v(oo). 

(ii) Using the relation between the mass-loss rates and the terminal wind ve- 
locities fitted by Winters et al. (2003): 

log lo (M/M yr^ 1 ) = -7.40 + 1 log^M/kms" 1 ). (1) 

With ascent of star along giant branch, the mass-loss rate becomes higher 
and higher, which results in the higher -y(oo). However, Eq. (1) is valid for M 
close to 10~ 6 M Q yr" 1 . For a mass-loss rate higher than 1O _6 M yr" 1 , Eq. (1) 
gives too high v(oo). Based on the models of Winters et al. (2000), we assume 
v(oo) = 30km s -1 if v (00) > 30km s" 1 . 

AM c ^ t D : The critical ignition mass of the nova depends mainly on the mass of 
accreting WD, its temperature and material accreting rate. Following Lii et 
al. (2006), we use AM™s by Eq. (Al) of Nelson et al. (2004) and Eq. (16) of 
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Table 1 

Parameters of the models for the population of SSs and SSSs in SySs. The first 
column gives the model number. Columns 2, 3 and 4 show the algorithm of the 
common envelope, the terminal velocity and the critical ignition mass, respectively. 

Eq.(Al) in Nelson et al. (2004) 
Eq.(Al) in Nelson et al. (2004) 
Eq.(Al) in Nelson et al. (2004) 
Eq.(16) in Yungelson et al. (1995) 



Yungelson et al. (1995). As far as whole populations of SySs are considered, 
the former gives lower AM*f than the latter (Lii et al., 2006). 



2.2 X-ray sources 

van den Heuvel et al. (1992) explained super-soft X-ray emissions by steady 
nuclear burning of hydrogen accreted onto WD. Most of the known SySs are 
sufficiently bright X-ray sources. As mentioned in Introduction, Miirset et 
al. (1997) divided X-ray sources into SSSs, SSs and relatively hard X-ray 
resources. It is well known that the SSSs in SySs are hot WDs whose photo- 
spheres can produce sufficiently hard ihv « 0.2kev) photons. For SSs in SySs, 
Jordan et al. (1994), Formiggini et al. (1995) and Miirset et al. (1997) sug- 
gested that they originate from the collision of two stellar winds. During some 
strong thermonuclear runaways, WD can eject some materials with high ve- 
locity (~ 1000 km s _1 ). Cool giant usually has a high mass-loss rate with a low 
velocity (5 — 30 km s _1 ). Therefore, a violent collision of the stellar winds is 
unavoidable in the strong thermonuclear runaways of SySs. An eruption of the 
recurrent nova RS Oph on February 12, 2006 provided the opportunity to per- 
form comprehensive X-ray observations. Nelson et al. (2008) showed its X-ray 
spectroscopy from 0.33 to 10 kev. The spectra indicated a collisionally dom- 
inated plasma with a broad range of temperature and an energy-dependent 
velocity structure. Sokolooski et al. (2006) and Orlando et al. (2009) suggested 
that most of the early X-ray emission of the 2006 outburst of RS Oph origi- 
nates from the interaction between the high-velocity ejecta from WD and the 
circumstellar medium which results mainly from the stellar wind of the cool 
giant in SyS. 

The high-velocity ejecta is a key factor for producing soft X-ray emissions in 
SySs. According to Yaron et al. (2005), Lii et al. (2006) divided the thermonu- 
clear runaways in SySs into two varieties: weak symbiotic nova in which most 
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of the accreted mater is deposited at the surface of the WD accretors; strong 
symbiotic nova in which WD accretors eject the majority of the accreted mater 
via high velocity wind. We assume SSs in SySs stem from the strong symbiotic 
novae, and the steady hydrogen burning 'ordinary' SySs and weak symbiotic 
novae result in SSSs. 

The strength of a thermonuclear runaway depends on the mass and the mass 
accretion rate of the WD. Using the results in Yaron et al. (2005), Lii et al. 
(2006) roughly defined the boundary between strong and weak symbiotic nova 
via the mass accretion rate and the mass of WD. We adopt the descriptions 
in Lii et al. (2006). Following Yungelson et al. (1996), we take the time which 
takes the WD to decline by 3 mag in its bolometric luminosity, t^ Q \, as the 
lifetime of SSSs. For SSs which originate from the interaction between the high- 
velocity ejecta from WD and the circumstellar medium, we take the duration of 
mass loss during strong symbiotic nova, t m \, as the their lifetime. One should 
note that £31,01 and £ m i are only a zero-order approximation. By a bilinear 
interpolation (Press et al. , 1992) of Table 3 in Yaron et al. (2005), t3 DO i and 
t m \ are calculated from the models in which the temperature of WD equals 
10 x 10 6 K. If the mass or the mass accretion rate of the WD in SSs are not 
in the range of the bilinear interpolation, they are taken as the most vicinal 
those. 



2. 3 Basic parameters of the Monte Carlo simulation 

We carry out binary population synthesis via Monte Carlo simulation tech- 
nique in order to obtain the properties of SSs and SSSs' population in SySs. 
For the population synthesis of binary stars, the main input model parameters 
are: (i) the initial mass function (IMF) of the primaries; (ii) the mass-ratio 
distribution of the binaries; (iii) the distribution of orbital separations; (iv) 
the eccentricity distribution; (v) the metallicity Z of the binary systems. 

We take IMF in Kroupa et al. (1993) as the primary mass distribution. The 
lower and the upper mass cut-offs of the primaries are O.8M and 8M Q , re- 
spectively. For the mass-ratio distribution of binary systems, we consider only 
a constant distribution(Kraicheva et al. , 1979; Mazeh et al., 1992; Goldberg 
& Mazeh, 1994). 

The distribution of separations is given by 

loga = 5X + l, (2) 

where X is a random variable uniformly distributed in the range [0,1] and a 
is in R Q . 



5 



In our work, the metallicity Z=0.02 is adopted. We assume that all binaries 
have initially circular orbits, and we follow the evolution of both components 
by BSE code, including the effect of tides on binary evolution (Hurley et al., 
2002). 

We assume that one binary with primary mass more massive than 0.8 M Q is 
formed annually in the Galaxy to calculate the birthrate of SSs and SSSs in 
SySs (Yungelson et al., 1994; Han et al., 1995a,b). 



3 Results 

We construct a set of models in which we vary different input parameters 
relevant to the population of SSSs and SSs in SySs. Table 1 gives all cases 
considered in this work. Many observational evidences showed that the ter- 
minal velocity of stellar wind v(oc) increases when a star ascends along the 
AGB (Olofsson et al, 2002; Winters et al, 2003; Bergeat & Chevallier, 2005). 
In this work, we take v(oc) calculated by Eq. (1) as the standard terminal 
velocity of stellar wind. 

We take 1 x 10 6 initial binary systems for each case. For every simulation with 
1 x 10 6 binaries, the relative errors for the symbiotic systems are lower than 
1%. Thus, 1 x 10 6 initial binaries appear to be an acceptable sample for our 
study. The main results of our study are given in Table 2 and in Figures 1-4. 

3.1 Population of X-ray sources 

According to Yungelson et al. (1995) and Lii et al. (2006), all progenitors of 
SySs pass through one of the three routes: 

(i) Channel I: unstable Roche lobe overflow (RLOF) with formation of a com- 
mon envelope; 

(ii) Channel II: stable RLOF; 

(hi) Channel III: formation of a white dwarf+red giant pair without RLOF. 
Table 2 gives the results of SSs' and SSSs' populations in SySs which have 
undergone channels I, II or III, respectively. 

In our model, SSs in SySs result from strong symbiotic novae. As Table 2 
shows, the Galactic occurrence rate of SSs in SySs is from about 2 (case 4) to 
20 (case 3) yr -1 . Channels I, II and III can produce SSs. According to Yaron 
et al. (2005), the typical duration of mass loss (t m i) for a strong outburst is 
about dozens of days. Therefore, SSs are transitory X-ray sources. SSSs in 
SySs result from weak symbiotic novae or steady hydrogen burning 'ordinary' 
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Table 2 

Different models of the population of SSSs and SSs in SySs. The first column gives 
case number according to Table 1. Columns 2, 3 and 4 give the occurrence rate of 
SSs which have undergone channel I, II and III, respectively. Column 5 shows the 
total occurrence rate of SSs. Columns 6, 7 and 8 give the numbers of SSSs which 
have undergone channel I, II and III, respectively. Column 9 shows the total number 
of SSSs. In these columns, the numbers in parentheses mean SSSs which originate 
from weak symbiotic novae. Columns 10, 11 and 12 show the occurrence rate of SSSs 
which originate from weak symbiotic novae and have undergone channel I, II and 
III, respectively. Column 13 shows the total occurrence rate of SSSs which originate 
from weak symbiotic novae. 
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SySs. For SSSs from the former, their number is between 46 (case 4) and 620 
(case 3), and their occurrence rate is from 2 (case 4) to 17 yr" 1 (case 2). For 
SSSs from the latter, their formation rate [between 0.01 (cases 1, 3 and 4) 
and 0.03 (case 2) yr -1 ] is very low and is not shown in Table 2. However, the 
lifetime of SSSs from steady hydrogen burning is very long, which results in 
large numbers [ between 2244 (case 4) and 5537 (case 2)]. Due to 7-algorithm 
which gives wide binary separation after common envelope evolution, SSSs 
in case 2 can be produced by channel I. While, SSSs in other cases do not 
undergo channel I. 

As the above descriptions, the input physical parameters (the algorithm of the 
common envelope, terminal velocity of the stellar wind v(oo) and the critical 
ignition mass AM^J 3 ) have great effects on the population of SSs and SSSs in 
SySs. For SSs, v(oo) is a key factor, and the Galactic occurrence rate of SSs in 
case 3 are about 5 times of those in case 1. For SSSs from weak symbiotic novae, 
AM^f affects their Galactic number and occurrence rate within a factor of 4; 
for SSSs from 'ordinary' SySs, the algorithm of the common envelope affects 
their Galactic number and formation rate within a factor of 2.4. A detailed 
analysis for the effects of these parameters can be seen in Lii et al. (2006). 

According to Lii et al. (2006), the Galactic numbers of SySs in the cases with 
the same input physical parameters are 5400 in case 1 and 4300 in case 4. The 
number ratios of the X-ray sources including SSs and SSSs to SySs are ~ 50% 
in case 1 and ~ 60% in case 4. v(oo) in case 3 is unfavorable for the formation 
of stable hydrogen burning 'ordinary' SySs (Lii et al., 2006), the number ratio 
is about 20%. Miirset et al. (1997) found that 60% of the known Galactic SySs 
are sufficiently X-ray bright to be detected. This result is agreement with our 
result. On comparing our results with observations, one should note that the 
number of X-ray sources depends greatly on the lifetime of the thermonuclear 
runaways as X-ray sources. 

In addition, it is very difficult to detect SSs and SSSs. The main reason is that 
dominant interaction with interstellar matter for X-rays between 0.1 keV and 
10 keV is photoelectric absorption, and the cross section increases rapidly at 
lower energies, scaling approximately as E~ 3 (McCammon & Wilton, 1990). 
According to Motch et al. (1994), SSSs have to be closer than ~ 2 kpc to 
the Sun in order to be detected. This means that only about ^5 of SSSs 
can be observed. There are 4 SSSs whose distances from the Sun are shorter 
than 2.6 kpc in Miirset et al. (1997). Hence, the Galactic number of SSSs 
in SySs is about 200 if 4 SSSs in SySs are a completed sample within 2.6 
kpc from the Sun. However, the orbital motion, the magnetic field and the 
gravitational influence of the companions in SySs can complicate the structure 
of circumstellar material. The nova ejecta and winds from the cool components 
in SySs act as shielding agents. Therefore, 4 SSSs are a uncompleted sample. 
A detailed work on X-ray spectrum absorptions by circumstellar material and 
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6 2 

lo § • p orb( da y s ) 

Fig. 1. Number distributions of X-ray sources in the Galactic SySs as a function of 
orbital periods (all numbers are normalized to 1). The solid, dashed and dotted lines 
represent SSs, SSSs from weak symbiotic novae and from steady hydrogen burning 
'ordinary' SySs, respectively. 

interstellar medium is very difficult and beyond the scope of this work. 



3.2 Properties of X-ray sources 



In this subsection, we describe potentially observable physical quantities of 
SSSs and SSs in SySs. 

Figure 1 gives the distributions of X-ray sources in the Galactic SySs over 
orbital periods. As we can see from case 1, SSSs have periods between 100 
and 100000 days. SSSs from steady hydrogen burning 'ordinary' SySs and 
from weak symbiotic novae have strong peaks at 7000 days and 11000 days, 
respectively. While SSs have orbital periods from 10 to 1000000 days, with 
a typical peak at 36000 days. According to Yungelson et ah (1995) and Lii 
et al. (2006), SySs via channel I have short orbital periods, while they via 
channel III have long orbital periods. As Table 2 shows, SSs can be produced 
via channels I, II and III, but SSSs can be done via channels II and III (except 
case 2). Therefore, SSs span wider range of orbital periods than SSSs. 

Figure 2 shows the distributions of numbers of the modeled X-ray sources 
in SySs each as a function of WD mass. The distribution for case 2 is wider 
than that in other cases. The main reason is that 7 algorithm for common 
envelope can form SySs with low-mass WD accretors (Lii et ah, 2006). For 
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Fig. 2. Number distributions of X-ray sources in the Galactic SySs as a function of 
WD mass (all numbers are normalized to 1). The solid, dashed and dotted lines rep- 
resent SSs, SSSs from weak symbiotic novae and steady hydrogen burning 'ordinary' 
SySs, respectively. 
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Fig. 3. Gray-scale maps of WD's mass vs. orbital period distributions for SSs and 
SSSs in SySs for case 1. The gradations of gray-scale correspond to the regions 
where the number density of systems is, respectively, within 1 - 1/2, 1/2 - 1/4, 1/4 
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all models, the distributions have typical masses in the O.5-O.8M range, but 
the peaks of the distributions are different from each type of X-ray sources. 
For example, in case 1, SSs have a peak at about O.7M . While, SSSs have 
a peak at about O.6M , which is in agreement with the results in Yungelson 
et al. (1996). According to Lii et al. (2006) and Yaron et al. (2005), massive 
WD favors to produce strong thermonuclear outburst. Mikolajewska (2003) 
gave the distribution of the measured masses of the WD in SySs. The peak 
is at about 0.5 M & regardless of Mikolajewska (2003) estimating low limits 
(Lii et al., 2006). Iben (2003) showed average mass of WD in symbiotic novae 
systems which are SSs is 0.73M Q . It looks likely that WD accretors in SSs have 
masses slightly larger than those in SSSs. Our results are compatible with the 
observations. 



Figure 3 shows the distributions of SSs and SSSs in case 1, in the " WD's 
mass - orbital period" plane. The ranges of orbital periods and WD's masses 
in SSs are wider than those in SSSs. The main reason is that producing SSSs 
needs higher mass accretion rates and the critical ignition mass AM^J 3 than 
producing SSs. As the top panel in Figure 3 shows, the distribution in SSs is 
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Fig. 4. Number distributions of X-ray sources in the Galactic SySs as a function of 
the mass accretion rate of WD (all numbers are normalized to 1). The solid, dashed 
and dotted lines represent SSs, SSSs from weak symbiotic novae and from steady 
hydrogen burning 'ordinary' SySs, respectively. 

cut into two regions. The low region represents these SSs which have undergone 
channel I, while the top region is for SSs which have undergone channels II 
and III. 



The distributions of numbers of X-ray sources in the Galactic SySs over the 
mass accretion rate of WD, M acc , is given in Figure 4. According to our model, 
M acc determines the strength of the hydrogen burning. For low M acc , the 
accreting WD experiences the strong thermonuclear outburst which results in 
SSs. For high M acc , the accreting WD can undergo stable hydrogen burning. 
Therefore, the peaks of magnitude of M acc in Figure 4 are 10~ 9 , 10~ 8 and 
10~ 7 Af Q yr~ 1 for SSs, SSSs from weak symbiotic novae and from 'ordinary' 
SySs, respectively. In Iben (2003), theoretically estimated M acc of five SSs is 
between 10~ 9 and 10~ 8 M Q yr~ 1 , which is in good agreement with our results. 



4 Conclusions 



Assuming that SSs in SySs result from the violent collisions of the stellar 
winds during strong symbiotic novae, and SSSs from weak symbiotic novae or 
steady hydrogen burning SySs, we investigate the Galactic population of SSs 
and SSSs in SySs. The Galactic occurrence rates of SSs and SSSs are from ~ 
2 to 20 yr _1 , and ~ 2 to 17 yr" 1 , respectively. The numbers of X-ray sources 
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in SySs range from 2390 to 6120. The percentage of the X-ray sources in all 
SySs is between ~ 20% and 60%, which is in agreement with the observations. 

In the present paper we do not model the X-ray spectra emissions of SSs 
and SSSs, and do not consider the interstellar medium's absorptions for X- 
ray spectra. Our work only shows some primary results. It is necessary for 
a detailed model of SSs and SSSs to simulate X-ray spectra emissions and 
interstellar medium's absorptions. We will consider them in the next paper. 
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